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1. Introduction 

In 1998, at the QCD workshop in Paris, Rolf Baier asked me whether jets could be mea- 
sured in Au+Au collisions because he had a prediction of a QCD medium-effect (energy loss via 
soft gluon radiation induced by multiple scattering [[I]]) on color-charged partons traversing a hot- 
dense-medium composed of screened color-charges [Qj. I told him [Q] that there was a general 
consensus [Q] that for Au+Au central collisions at y/SNN = 200 GeV, leading particles are the only 
way to find jets, because in one unit of the nominal jet-finding cone, Ar = y 7 (At]) 2 + (A0) 2 , there 
is an estimated 7iAr 2 x j^'-jjj- ~ 375 GeV of energy !(!) The good news was that hard-scattering 
in p-p collisions was originally observed by the method of leading particles and that these tech- 
niques could be used to study hard-scattering and jets in Au+Au collisions. In fact, in several 
recent talks [||, ^, ^ and in talks from earlier years [^, |fl, some as long ago as 1979 I have been 
on record describing "How everything you want to know about jets can be found using 2-particle 
correlations". This past year, I had to soften the statement to almost everything because we found 



by explicit calculation in PHENIX [10] that the two-particle opposite-side correlation is actually 
quite insensitive to the fragmentation function — overturning a belief dating from the seminal paper 
of Feynman, Field and Fox in 1977 [pl|]. However, we also found that the opposite-side correlation 
function is sensitive to the ratio of the transverse momentum of the away-side jet (prj to that of 
the trigger-side jet {pj t ) and thus provides a way to measure the relative energy loss of the two jets 
from a hard-scattering which escape from the medium in an A+A collision. 

2. Status of theory and experiment, circa 1982 

Hard-scattering was visible both at the ISR and at FNAL fixed-target-energies via inclu- 
sive single particle production at large pj > 2-3 GeV/c [|l2|]. Scaling and dimensional argu- 



ments [|13J, [14j, |15|, |16p for plotting data revealed the systematics and underlying physics. The 
theorists had the basic underlying physics correct; but many (inconvenient) details remained to be 
worked out, several by experiment. The transverse momentum imbalance of outgoing parton-pairs, 



the "kj -effect", was discovered by experiment [17, 18], and clarified by Feynman and collabora 



tors [11]. The first modern QCD calculation and prediction for high pj single particle inclusive 
cross sections, including non-scaling and initial state radiation was done in 1978, by Jeff Owens 
and collaborators [ |l9| ] under the assumption that high pj particles are produced from states with 
two roughly back-to-back jets which are the result of scattering of constituents of the nucleons 
(partons). 



The overall p-p hard-scattering cross section in "leading logarithm" pQCD [gOp is the sum 
over parton reactions a + b — > c + d (e.g. g + q — > g + q) at parton-parton center-of-mass (cm.) 
energy \/S: 

<±q<2A"2acos U* dsdydcosU* s ~r 2x\X2 

where f a {x\), fb(x%), are parton distribution functions, the differential probabilities for partons a 
and b to carry momentum fractions x\ and X2 of their respective protons (e.g. ufa)), and where 6* 
is the scattering angle in the parton-parton cm. system. The parton-parton cm. energy squared 
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is S = X1X2S, where y/s is the cm. energy of the p-p collision. The parton-parton cm. system 
moves with rapidity y = I /21n(jci/x2) in the p-p cm. system. The quantities f a (x\) and fyfa), the 
"number" distributions of the constituents, are related (for the electrically charged quarks) to the 
structure functions measured in Deeply Inelastic lepton-hadron Scattering (DIS), e.g. 

F l (x,Q 2 ) = W e2 aU x >Q 2 ) and F2(x,Q 2 )=x J £e 2 a f a (x,Q 2 ) (2.2) 

a a 

where e a is the electric charge on a quark. 

The Mandelstam invariants s, t and u of the constituent scattering have a clear definition in 
terms of the scattering angle 6* in the constituent cm. system: 




(2.3) 



The transverse momentum of a scattered constituent is: 



Pt=P*t = 




(2.4) 



and the scattered constituents c and d in the outgoing parton-pair have equal and opposite momenta 
in the parton-parton (constituent) cm. system. A naive experimentalist would think of Q 2 = —t for 
a scattering subprocess and Q 2 = s for a Compton or annihilation subprocess. 

Equation ^j] gives the pj spectrum of outgoing parton c, which then fragments into a jet 
of hadrons, including e.g. k°. The fragmentation function (z) is the probability for a k° to 
carry a fraction z, = p 71 / p c of the momentum of outgoing parton c. Equation ^j] must be summed 
over all subprocesses leading to a 7T° in the final state weighted by their respective fragmentation 
functions. In this formulation, f a (xi), fb{xz) and Df(z) represent the "long-distance phenom- 
ena" to be determined by experiment; while the characteristic subprocess angular distributions, 
£ a *(cos0*) (see Fig. [I]) and the coupling constant, a s (Q 2 ) = 25 ^(qi/a 1 ) ' 316 fundamental predic- 
tions of QCD [21,^2|] for the short-distance, large-<2 2 , phenomena. When higher order effects are 
taken into account, it is necessary to specify factorization scales /I for the distribution and frag- 
mentation functions in addition to renormalization scale A which governs the running of a s {Q 2 ). 
As noted above, the momentum scale Q 2 ~ p\ for the scattering subprocess, while Q 2 ~ s for a 
Compton or annihilation subprocess, but the exact meaning of Q 2 and \i 2 tend to be treated as 
parameters rather than as dynamical quantities. 

Due to the fact (which was unknown in the 1970's) that jets in 4n calorimeters at ISR ener- 
gies or lower are invisible below Vj ~ Ej < 25 GeV [25], there were many false claims of jet 
observation in the period 1977-1982. This led to skepticism about jets in hadron collisions, par- 
ticularly in the USA [pq|. A 'phase change' in belief-in-jets was produced by one UA2 event at 
the 1982 ICHEP in Paris [ |23| ] (Fig. ||), which, together with the first direct measurement of the 
QCD constituent-scattering angular distribution, L ab (cosd*) (Eq. |Q| ), using two-particle correla- 
tions [24], presented at the same meeting (Fig. II]), gave universal credibility to the pQCD descrip- 



tion of high pj hadron physics [20, 27 



4 



Review of hard scattering and jet analysis 



Michael J. Tannenbaum 
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Figure 1: a) (left 3 panels) CCOR measurement [23, 24 1 of polar angular distributions of n° pairs with net 
Pt < 1 GeV/c at mid-rapidity in p-p collisions with ^fs = 62.4 GeV for 3 different values of %n invariant 
mass M KK . b) (rightmost panel) QCD predictions for E^fcos 0*) for the elastic scattering of gg, qg, qq', qq, 
and qq with OC s (Q 2 ) evolution. 





(a) 



(b) 



Figure 2: UA2 jet event from 1982 ICHEP [|23j]. a) event shown in geometry of detector, b) "Lego" plot 
of energy in calorimeter cell as a function of angular position of cell in polar (6) and azimuthal (<t>) angle 
space. 



3. Mid-rapidity pj spectra from p-p collisions — *r-scaling 



Equation 2.1 leads to a general l xj -scaling' form for the invariant cross section of high-/?r 
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particle production: 



J 3 Q 



1 



n(x T ,J~s) 

Pt 



F(x T ) 



1 



-G{x T ) , 



(3.1) 



where xj = 2pj j yfs. The cross section has two factors, a function F{xj) (G(xt)) which 'scales', 
i.e. depends only on the ratio of momenta, and a dimensioned factor, l/p^ XT '^ (1/ y/s^ 7 '^), 
where n{xj, yfs) equals 4 in lowest-order (LO) calculations, analogous to the l/q 4 form of Ruther- 
ford Scattering in QED. The structure and fragmentation functions are all in the F{xt) (G(xt)) 
term. Due to higher-order effects such as the running of the coupling constant, a s (Q 2 ), the evo- 
lution of the structure and fragmentation functions, and the initial-state transverse momentum kj, 
n(xr,y/s) is not a constant but is a function of x T , y/s. Measured values of n(x T ,y/s) for 71° in p-p 
collisions are between 5 and 8 [Q]. 

The scaling and power-law behavior of hard scattering are evident from the yfs dependence 
of the pt dependence of the p-p invariant cross sections. This is shown for nonidentified charged 
hadrons, (h + + h~)/2, in Fig. ||a. At low pj < 1 GeV/c the cross sections exhibit a "thermal" 
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Figure 3: a) (left) Ed 3 a{pT)/d i p at mid-rapidity as a function of y/s in p-p collisions j2^]. b) (right) 
log-log plot of V?(GeV) 63 x Ed 3 a/d 3 p vs x r = 2p T /yfs [M. 



exp (—6pr) dependence, which is largely independent of y/s, while at high pj there is a power-law 
tail, due to hard scattering, which depends strongly on yfs. The characteristic variation with y/s at 



high pj is produced by the fundamental power-law and scaling dependence of Eqs. |2JJ, |3.1[ This 
is best illustrated by a plot of 

d 3 a 



n(XT,y/s) 



XE 



d 3 p 



G(x T ) 



(3.2) 



as a function of xj, with n(xr,y/s) = 6.3, which is valid for the xj range of the present RHIC 
measurements (Fig. |3]b). The data show an asymptotic power law with increasing xj. Data at 
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a given y/s fall below the asymptote at successively lower values of xj with increasing yfs, cor- 
responding to the transition region from hard to soft physics in the pj region of about 2 GeV/c. 
Although xj -scaling provides a rather general test of the validity QCD without reference to details, 
the agreement of the PHENIX measurement of the invariant cross section for 71° production in p-p 



collisions at yfs = 200 GeV Q with NLO pQCD predictions over the range 2.0 < p T < 15 GeV/c 
(Fig. fib is, nevertheless, impressive. 
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Figure 4: (left) PHENIX O 



n° invariant cross section at mid-rapidity from p-p collisions at y/s = 200 GeV, 
together with NLO pQCD predictions from Vogelsang [31, 32[] . a) The invariant differential cross section for 
inclusive n° production (points) and the results from NLO pQCD calculations with equal renormalization 
and factorization scales of pr using the "Kniehl-Kramer-Potter" (solid line) and "Kretzer" (dashed line) sets 
of fragmentation functions, b) The relative statistical (points) and point-to-point systematic (band) errors. 
c,d) The relative difference between the data and the theory using KKP (c) and Kretzer (d) fragmentation 
functions with scales of pj/2 (lower curve), pj, and 2pj (upper curve). In all figures, the normalization 
error of 9.6% is not shown, (right) e) p-p data from a) multiplied by the nuclear thickness function, T^a, 
for Au+Au central (0-10%) collisions plotted on a log-log scale (open circles) together with the measured 
semi-inclusive 7T° invariant yield in Au+Au central collisions at y/Im = 200 GeV [p3|| 



3.1 The importance of the power law 

A log-log plot of the 71° spectrum from Fig. |]a in p-p collisions, shown in Fig. ||e along with 
corresponding data from Au+Au collisions ||33|], illustrates that the inclusive single particle hard- 
scattering cross section is a pure power law for pj > 3 GeV/c. The invariant cross section for 71° 
production can be fit to the form 

Ed 3 a /dp 3 oc p T n (3.3) 
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with n = 8.10 ±0.05 [33]. It is important to emphasize that this n is different from the n(xT,\fs) 
in Eqs. [O] , |3.2| . The power n in Eq. [53] measures the pure power law shape of the cross section 



at a fixed y/s represented by the function G(xj) in Eqs. "SA , 3j , while n(xr, \fs) represents the pj 
dependence at fixed x T when y 7 ! is varied. Clearly, from Eq. 3T and Fig. ^p, the simple power n 



(Eq. p3[ ) is greater than or equal to n(xr, \/s) (Eq. 3.2). 

The steeply-falling power-law spectrum at a given y 7 ? has many important consequences for 
single particle inclusive measurements of hard-scattering. The most famous properties in this re- 



gard are the "Bjorken parent-child relationship" [34] and the "leading-particle effect", which also 
goes by the unfortunate name "trigger bias" [35, [34]. These will be discussed below in the section 
on correlations. The power-law also makes the calculation of the inclusive photon spectrum from 
the decay 71° — > y+ y very easy, but nevertheless very precise [36, [57|], when expressed as the ratio 
of photons from 71° to 7T° at the same pj : 



7T ,; 



(p T )=2/(n-\) 



(3.4) 



Similarly, the inclusive electron spectrum from internal or external conversion of these photons has 
a simple formula when expressed as the ratio to 7C° at the same pj : 



7T' 



(pt) 



(e- + e+) 



2tt u 



(p 



82 



T = 



7^0 



(3.5) 



where 82/2 = Dalitz (internal conversion) branching ratio per photon and t/Xo is the thickness of 



the external converter in radiation lengths (Xq) [iO, 41] 



4. Measurement of the medium effect in A+A collisions with hard-scattering by 
comparison to baseline measurements in p-p and d+A collisions 

Since hard scattering is point-like, with distance scale l/pr < 0.1 fm, the cross section in 
p+A (B+A) collisions, compared to p-p, should be simply proportional to the relative number of 



possible point-like encounters Q38J], a factor of A (BA) for p+A (B+A) minimum bias collisions. 
For semi-inclusive reactions in centrality class / at impact parameter b, the scaling is proportional 
to TAB{b), the overlap integral of the nuclear thickness functions [39], where (Tab) / averaged over 
the centrality class is: 



(T, 



AB/f 



f f T AB (b) 



d 2 b 



= (Ncqu)/ 

1 _ e -<*NN T AB (b) \ Onn 



(4.1) 



and where (N co ii ) / is the average number of binary nucleon-nucleon inelastic collisions, with cross 
section Onn, in the centrality class /. This leads to the description of the scaling for point-like pro- 
cesses as binary-collision (or N co ii) scaling. This description is convenient, but confusing, because 
the scaling has nothing to do with the inelastic hadronic collision probability, it is proportional only 
to the geometrical factor (Tab) / (Eq. |Q| ). 

Effects of the nuclear medium, either in the initial or final state, may modify the point-like 
scaling. This is shown rather dramatically in Fig. ||e where the Au+Au data are suppressed relative 
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to the scaled p-p data by a factor of ~ 4 — 5 for pj > 3 GeV/c. A quantitative evaluation of the 
suppression is made using the "nuclear modification factor", Rab, the ratio of the measured semi- 
inclusive yield to the point-like scaled p-p cross section: 

Rab= dN p AB _ dN p AB 

(T AB ) f x dofi N (N co „) f x dNfi N 

where dN AB is the differential yield of a point-like process P in an A + B collision and da^ N is the 
cross section of P in an AW (usually p-p) collision. For point-like scaling, Rab = 1- 

While the suppression of n° at a given pj in Au+Au compared to the scaled p-p spectrum 
may be imagined as a loss of these particles due to, for instance, the stopping or absorption of a 
certain fraction of the parent partons in an opaque medium, it is evident from Fig. ||e that an equally 
valid quantitative representation can be given by a downshift of the scaled p-p spectrum due to, for 
instance, the energy loss of the parent partons in the medium — a particle with p'j in the scaled p-p 
spectrum is shifted in energy by an amount S(pr) to a measured value pj = p'j — S(pj) in the 
Au+Au spectrum [43]. The fact that the Au+Au and reference p-p spectra are parallel on Fig. [|e 
provides graphical evidence that the fractional pj shift in the spectrum, S(pj)/ Pt is a constant for 
pj > 3 GeV/c, which, due to the power law, results in a constant ratio of the 71° spectra Raa(pt) 
as shown in Fig. |5[ 
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Figure 5: PHENIX measurement of nuclear modification factor Raa for identified n° and-non identified 
charged hadrons (h + +h~)/2 for central (0-10%) Au+Au collisions at ^/sim = 200 GeV Q. 

The nuclear modification factors are clearly different for n° and (h + + h~)/2 for pj < 6 GeV/c 
in Fig. |[ Is it possible to tell whether one or both of these reactions obey QCD? 

4.1 xj scaling in A+A collisions as a test of QCD 

If the production of high-^y particles in Au+Au collisions is the result of hard scattering ac- 
cording to pQCD, then xj scaling should work just as well in Au+Au collisions as in p-p collisions 
and should yield the same value of the exponent «(jtr,\A)- The only assumption required is that 



the structure and fragmentation functions in Au+Au collisions should scale, in which case Eq. 3.2 
still applies, albeit with a G(xj) appropriate for Au+Au. In Fig. ^, n(xj ,^/snn) in Au+Au is 
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derived from Eq. 3.2, for peripheral and central collisions, by taking the ratio of Ed 3 a/dp 3 at a 
given xj for Js^m = 130 and 200 GeV, in each case. The 7T°'s exhibit xj scaling, with the same 
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Figure 6: Power-law exponent n{xr) for 7T° and h spectra in central and peripheral Au+Au collisions at 
y/s^j = 130 and 200 GeV @]. 

value of n = 6.3 as in p-p collisions, for both Au+Au peripheral and central collisions, while the 
non-identified charged hadrons Aj-scale with n = 6.3 for peripheral collisions only. Notably, the 
(h + +h~)/2 in Au+Au central collisions exhibit a significantly larger value of n{xj, \fs), indicat- 
ing different physics, which will be discussed below. The xj scaling establishes that high-pr- 71 
production in peripheral and central Au+Au collisions and (h + + h~)/2 production in peripheral 
Au+Au collisions follow pQCD as in p-p collisions, with parton distributions and fragmentation 
functions that scale with xj, at least within the experimental sensitivity of the data. The fact that 
the fragmentation functions scale for 71° in Au+Au central collisions indicates that the effective 
energy loss must scale, i.e. S(pt)/pt = is constant, which is consistent with the parallel spectra 
on Fig. ||e and the constant value of Raa as noted in the discussion above. 

The deviation of (h + + hr)/2 from xj scaling in central Au+Au collisions is indicative of and 
consistent with the strong non-scaling modification of particle composition of identified charged- 
hadrons observed in Au+Au collisions compared to that of p-p collisions in the range 2.0 < pr < 
4.5 GeV/c, where particle production is the result of jet-fragmentation. As shown in Fig. [7|-(left) 
the p/n + and p/n~ ratios as a function of pj increase dramatically to values ~1 as a function 



of centrality in Au+Au collisions at RHIC Q45p which was totally unexpected and is still not fully 
understood. Interestingly, the p and p in this p T range appeal - to follow the N co u scaling expected 
for point-like processes (Fig ^-(right)), while the n° are suppressed, yet this effect is called the 
'baryon anomaly', possibly because of the non-^T- scaling. An elegant explanation of this effect as 
due to coalescence of quarks from a thermal distribution [[if], 47, ^8j, which would be prima facie 
evidence of a Quark Gluon Plasma, is not in agreement with the jet correlations observed in both 
same and away-side particles associated with both meson and baryon triggers [^] (see discussion 



of Fig. 24 below). 



4.2 Direct photon production 

Direct photon production is one of the best reactions to study QCD in hadron collisions, since 
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GeV [^5j. Open (filled) points are for IV s1 (7t°), respectively, (right) Invariant yield of p and p, from the 
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there is direct and unbiased access to one of the interacting constituents, the photon, which can be 



measured to high precision, and production is predominantly via a single subprocess [50]: 



g + q^J+q , (4.3) 

with q + q — ► J + g contributing on the order of 10%. However, the measurement is difficult ex- 
perimentally due to the huge background of photons from 71° — > 7 + 7 and rj — ► 7 + 7 decays. This 



background can be calculated using Eq. and can be further reduced by 'tagging' — eliminating 
direct-photon candidates which reconstruct to the invariant mass of a tt° when combined with 
other photons in the detector, and/or by an isolation cut — e.g. requirement of less than 10% ad- 
ditional energy within a cone of radius Ar = (At]) 2 + (A0) 2 = 0.5 around the candidate photon 
direction — since the direct photons emerge from the constituent reaction with no associated frag- 
ments. 

The exquisite segmentation of the PHENIX Electromagnetic calorimeter (At7 x A(j) ~ 0.01 x 
0.01) required in order to operate in the high multiplicity environment of RHI collisions also pro- 
vides excellent 7 and 71° separation out to pj ~ 25 GeV/c. This will be useful in making spin- 
asymmetry measurements of direct photons in polarized p-p collisions for determination of the 



gluon spin structure function pip, but, in the meanwhile, has provided a new direct photon mea- 
surement in p-p collisions which clarifies a longstanding puzzle between theory and experiment in 
this difficult measurement. In Fig. |8|-(left) the new measurement of the direct photon cross sec- 
tion in p-p collisions at ^f~s = 200 GeV from PHENIX []52|] is shown compared to a NLO pQCD 
calculation, with excellent agreement for pj > 3 GeV/c. This data has resolved a longstanding 



discrepancy in extracting the gluon structure function from previous direct photon data [ ]53| , |54j] 
(see Fig. ||-(right)) by its agreement with ISR data and the theory at low xj. 

4.3 xj -scaling in direct photon, jet and identified proton production in p-p collisions 

The new direct photon measurement also shows nice xj scaling with previous measurements 
(Fig. 0-(left)) with a value n(xr , y/s) = 5 .0. This is closer to the asymptotic value of n{xj , yfs) = 4 
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Figure 8: (left) (a) Direct photon spectra [|2|] with NLO pQCD calculations for three theory scales, ju. 
Brackets around data points show systematic errors, (b) Comparison to the NLO pQCD calculation for 
ji = Pt , with upper and lower curves for p, = pr/2 and 2pr- (right) Ratio of direct photon measurements 
to theoretical calculation [|3| |4j 



than the 71° measurements in this range of pr, y/s but is still not as close as the n{xr,y/s~) = 4.5 
from jet measurements [^] at the Tevatron (Fig. ||-(left)). 

A new measurement by STAR [56] of xj scaling of identified n ± , p and p in p-p collisions 
at y/s = 200 GeV in comparison to previous measurements (Fig. pi]) gives n(xT,y/s) = 6.8 ±0.5 
for pions in agreement with the PHENIX measurement [Q] and provides the first measurement of 
xj scaling of p and p in p-p collisions, with n(xr, y/s) = 6.5 ± 1.0 in agreement with the value for 
pions. This result shows that p and p are produced by fragmentation of hard-scattered partons in 
p-p collisions for pj > 2 GeV/c, which contradicts a recent proposal [57] to explain the 'baryon 
anomaly' in A+A collisions as due to the possibility that protons and pions in the range 2.0 < 
Pt < 4.5 GeV/c in p-p collisions are produced by different mechanisms. As shown on the inset in 
Fig. 1C , the pion and proton spectra follow transverse mass scaling for mj < 2 GeV/c 2 in both p-p 
and d+Au collisions, suggesting the transition region from soft to hard process domination occurs 
at pr ~ 2 GeV/c in these collision systems. 



4.4 The state of jet-suppression measurements in Au+Au and d+Au collisions at RHIC 

The state of Raa measurements at RHIC is beautifully summarized in Fig. |TI] where the nuclear 
modification factor is the same for n° and tj in Au+Au central collisions at ^/snn = 200 GeV, both 
are suppressed relative to point-like scaled p-p data by a factor of ~ 5 which appears to be constant 
for pj > 4 GeV/c, while the direct photons are not suppressed at all. Since the direct photons do not 
interact (strongly) with the medium, while the 71° and v\ are fragments of outgoing hard-scattered 
partons which do interact with the medium, this plot proves that the suppression is a medium effect. 
The curve on the plot shows a theoretical prediction [59] from a model of parton energy loss. The 
model assumes an inital parton density dN/dy = 1200, which corresponds to an energy density 
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Figure 9: (left) xj scaling J52| ] of direct photon data in p-p and p-p collisions. The quantity plotted is 
(\/s) n x Ed^a /dp i {xj) with n = 5.0. (right) xj scaling of jet cross sections measured in p-p collisions by 
CDF and DO J55]]. The quantity plotted is the ratio of p\ times the invariant cross section as a function of 
xt for y/s = 630 and 1800 GeV. Note that the theory curves are plotted in the same way in order to avoid as 
much as possible uncertainties from the various parton distribution functions used. 



of approximately 15 GeV/fm 3 . The theory curve appears to show a reduction in suppression with 
increasing pj, while, as noted above, the data appear to be flat to within the errors, which clearly 
could still be improved. 

It is unreasonable to believe that the properties of the medium have been determined by a 
theorist's line through the data which constrains a few parameters of a model. The model and 
the properties of the medium must be able to be verified by more detailed and differential mea- 
surements. All models of medium induced energy loss [ ]60| ] predict a characteristic dependence of 
the average energy loss on the length of the medium traversed. This is folded into the theoretical 
calculations with added complications that the medium expands during the time of the collision, 
etc [61]. In an attempt to separate the effects of the density of the medium and the path length 
traversed, PHENIX [^, |62|] has studied the dependence of the 71° yield as a function of the an- 
gle (Aip) to the reaction plane in Au+Au collisions (see Fig. 12). For a given centrality, variation 
of A(p gives a variation of the path-length traversed for fixed initial conditions, while varying the 
centrality allows the initial conditions to vary. Clearly these data reveal much more activity than 
the reaction-plane-integrated Raa (Fig- 0) an( ^ mer it further study by both experimentalists and 
theorists. 

The point-like scaling of direct photon production in Au+Au collisions indicated by the ab- 
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Figure 10: x^-scaling of pions, protons and anti-protons as presented by STAR [ pq ] in comparison to 
previous measurements. The inset shows the m^-scaling of the invariant yields of charged pions and protons+ 
anti-protons in p-p and d+Au collisions. 
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Figure 11: Nuclear modification factor, Raa for direct photons, 7T° and tj in Au+Au central collisions at 
V^wv = 200 GeV 
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Figure 12: Nuclear modification factor, Raa(Pt)\a$, of 7T°'s as a function of pj for several values of cen- 
trality and angle to the reaction plane (A0) in Au+Au collisions at ^Aaw = 200 GeV [ j33[ |62[ ], The shaded 
region around R^a = 1 indicates the systematic error due to the reaction plane resolution correction. The 
arrows sketched on the ellipse indicate the different path-lengths traversed for centrality ~ 30 — 40% . 

sence of suppression in Fig. |TT] implies that there also should be no suppression of direct pho- 
tons in d+A collisions. The further implication is that the gluon structure function in a nucleus 
(^(x)) scales like A, i.e. R^(x) = g A (x)/Ag N \x) = 1, where g N (x) is the partem distribution func- 
tion of gluons in a nucleon, since it is known from measurements of deeply inelastic scattering 
of muons in nuclei that there is only a slight effect in the quark-structure functions, R%(x) = 
F£(x)/AF£(x)<\ @||, U (see Eq. However, until now there has been no direct measure- 



ment of the gluon structure function in nuclei. A first attempt in this direction has been presented 
by the PHENIX collaboration as a measurement of R c /a{pt) of direct photons in d+Au collisions 
at y/sjfij = 200 GeV (Fig. 0-(left)). 

At mid-rapidity there is a simple relationship between R^a for direct photon production and 



the structure function ratios [50] 



(4.4) 

The measurement is consistent with R A (x) = 1, but clearly the statistical errors must be improved by 



an order of magnitude before the data can be compared in detail to the theoretical prediction [pfl, 67] 
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Figure 13: (left)-Nuclear modification factor R c /a of direct photons (red solid circles) and neutral pions 
(open circles) measured in minimum bias d+Au collisions at ^/snn = 200 GeV [|65[]. (right) Theoretical 
prediction for gluon structure function ratio Rg(x) in Pb [ |6(| p7§ . Note that x ^ xj — /?r(GeV/c)/100 for 
comparing the left and right figures. 



(Fig. 13 -(right)) which is used in all calculations for RHI collisions. 



5. Correlations 



As noted above (section 3J), the steeply-falling power-law spectrum at a given y/s has many 
important and helpful consequences for single particle inclusive and two-particle correlation mea- 
surements of hard-scattering. The most famous properties in this regard are the "Bjorken parent- 
child relationship" [p4l] and the "leading-particle effect", which also goes by the unfortunate name 



"trigger bias" [B3, 34]. 



5.1 Why single particle inclusive measurements accurately measure hard-scattering — the 
leading-particle effect, also known as "trigger bias." 

Due to the steeply falling power-law transverse momentum (pj,) spectrum of the scattered 
partem, the inclusive single particle (e.g. it) pj t spectrum from jet fragmentation is dominated by 
fragments with large Zt, where Zt = Pt,/Pt, is the fragmentation variable. The joint probability for 
a fragment pion, with pj t = ZtpT,, originating from a partem with pj t = prjet is: 



d 2 On{PT,,Zt) 

p T ,dp T ,dzt 



PT,dp T , 

f q (pT,)xD«(z t ) 



(5.1) 



where f q (pT,) represents the final-state scattered-parton invariant spectrum do q / p^dpT, and D^(z t ) 
represents the fragmentation function. The first term in Eq. ^j] is the probability of finding a partem 
with transverse momentum pj t and the second term corresponds to the conditional probability that 
the parton fragments into a particle of momentum p T = z t pT,- A simple change of variables, p Tt = 
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PTt/zt, dpT,/dpT,\ Zt = 1/Zf , then gives the joint probability of a pion with transverse momentum pj t 
which is a fragment with momentum fraction Zt from a parton with pj t = PT,/zt- 

1 j =Jq{ ) xD q{Zt)X— ■ (5.2) 

p Tt dp Tl dz t z t q zf 

The pT t and zt dependences do not factorize. However, the pj t spectrum may be found by in- 
tegrating over all values of pj t from p Tt < p Tt < yfs/2, which corresponds to values of z.t from 
xj =2pr I 'y/sto 1: 

j =/ fq{ — ) D q{Zt)— ■ (5.3) 

PT,dp T , Jxt Zt q Zt 1 

Also, for any fixed value of p Tt one can evaluate the (z t {pT,)}> integrated over the parton spectrum: 

, , „ LUD",(Z,)UPT,/Z,)f 



Since the observed jt° spectrum is a power-law for 777; > 3 GeV/c, one can deduce from Eq. 53 



that the partonic pj t spectrum is also a power-law with the same power — this is the 'Bjorken parent- 
child relationship" [[34]]. If we take: 



do q 

PT,dp T , 

then 



f q (p Tt )=Apr" , (5.5) 



1 d ° % ~ fl AD«( Zt )( PT >^ dZt 



PT,dp Tl Jx T " q Zt zj 

■1 



l - f AD«(zt)zr 2 dz t , (5.6) 

'7* JXT 



Pt, -> x t 

where the last integral depends only weakly on pj t due to the small value of xj. Eq. ^6| also 
indicates that the effective fragmentation function for a detected inclusive single particle (with 
Pt,) is weighted upward in z.t by a factor z"~ 2 , where n is the simple power fall-off of the jet 



invariant cross section (i.e. not the n(xj,\fs) of Eq. 3.1). This is the so called "trigger bias" 
although it doesn't actually involve a hardware trigger. Any particle selected from an inclusive pj t 
spectrum will most likely carry a large fraction of its parent parton transverse momentum; and it 



was commonly accepted that this would define the hard scattering kinematics (s in Eq. 2.4 or pj t 



in Eq. 5.1 ) so that the jet from the other outgoing parton in the hard-scattered parton-pair would be 



unbiased [J73J, 11 1, so that its properties such as the fragmentation function and the fragmentation 



transverse momentum could be measured. 

5.1.1 Fragmentation Formalism — Single Inclusive 

For an exponential fragmentation function, 

D{z)=Be- bz , (5.7) 
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calculation of the "trigger bias" and the "parent-child" factor is straightforward [35]. The mean 
multiplicity of fragments in the jet is: 



( m )= [ D(z)dz= T (l-e- b ) 
Jo b 

and these fragments carry the total momentum of the jet: 

[ l zD(z)dz=^(l-e- b (\+b)) = \ 
Jo b z 



where the (z) per fragment is: 



The results are: 



B 



(m) 



(z) 



£zD{z)dz _ 1 
J l D(z)dz H 

b 2 

~ \ -e- b (l+b) ^ 
b(l-e- h ) 



l-e- b (l+b) 
l-e- b {l + b) 1 



b 



(5.8) 
(5.9) 

(5.10) 

(5.11) 
(5.12) 
(5.13) 



b{l-e- b ) b 

The mean multiplicity of charged particles in the jet is (m) « b, which is 8-10 at RHIC (see below). 



Substitution of Eq. £J\ into Eq. 5J5 for the pj t spectrum of the % gives: 

1 do n AB f l 



Pt, dp Tl p\ Jx Tt 



AB f l 

— / dztzT 2 exp -bz t 

Pt Jxt, 



which can be written as: 



~P L = wi~T inn ~ hbx Tl )-T{n -l,b)} 
Pt, dp T , Pt, b" 



where 



T(a,x) = 



t"- l e-'dt 



(5.14) 



(5.15) 



(5.16) 



is the Complementary or upper Incomplete Gamma function, and r(a,0) = T(a) is the Gamma 
function, where T(a) = (a — 1 ) ! for a an integer. 

A reasonable approximation for small xj values is obtained by taking the lower limit of 
Eq. 5.14 to zero and the upper limit to infinity, with the result that: 

1 da n T{n-\)AB 



Pt, dp Tt b n 1 p% 



(5.17) 



The parent-child ratio, the ratio of the number of % at a given pt, to the number of partons at the 
same p^ is just given by the ratio of Eq. |5.17| to Eq. at Pt, = Pt, '■ 



K 



Br(n-1) (m)r(n-l) 

\Pt,) = 

71° 



\n— 1 



W-2 



(5.18) 



18 



Review of hard scattering and jet analysis 



Michael J. Tannenbaum 



Similarly, the same substitutions inEq. 5.4 for (zt(pT,)) give: 

X£. 1 ex P -bzt 1 \T(n,bxT t ) - r(n,b)] 



(Zt{PT t )) 



Jl f dz t zr 2 exp -bz t b [T(n - 1 , ta r , ) - T(n - 1 , &)] b 



(5.19) 



This shows the "trigger-bias" quantitatively. The (zt{pT,)) of an inclusive single particle (e.g 
71°) with transverse momentum pj t , which is a fragment with momentum fraction z. t from a parent 
parton with pr t = pr,/zt (Eq 5.19 ), is n — 1 times larger than the unconditional (z) of fragmentation 



(Eq. |5.13| ) Q68Q. The prevailing opinion from the early 1970's until early this year was that although 
the inclusive single particle (e.g. pizero) spectrum from jet fragmentation is dominated by trigger 
fragments with large (z t ) ~ 0.7 — 0.8 the away-jets should be unbiased and would measure the 
fragmentation function, once the correction is made for (z t ) and the fact that the jets don't exactly 



balance pj due to the kj smearing effect [12] 



5.2 Almost everything you want to know about jets can be found using 2-particle 
correlations. 

The outgoing jet-pair of hard-scattering obeys the kinematics of elastic-scattering (of partons) 
in a parton-parton cm. frame which is longitudinally moving with rapidity y = 1 /21n(jci jxj) in the 
p-p cm. frame. Hence, the jet-pair formed from the scattered partons should be co-planar with the 
beam axis, with two jets of equal and opposite transverse momentum. Thus, the outgoing jet-pair 
should be back-to-back in azimuthal projection. It is not necessary to fully reconstruct the jets in 
order to measure their properties. In many cases two-particle correlations are sufficient to measure 
the desired properties, and in some cases, such as the measurement of the net transverse momentum 
of a jet-pair, may be superior, since the issue of the systematic error caused by missing some of 
the particles in the jet is not-relevant. Many ISR experiments provided excellent 2-particle corre- 
lation measurements [59]. However, the CCOR experiment [ f70| ] was the first to provide charged 
particle measurement with full and uniform acceptance over the entire azimuth, with pseudorapid- 
ity coverage —0.7 < tj < 0.7, so that the jet structure of high pj scattering could be easily seen 



and measured. In Fig. |14|a,b, the azimuthal distributions of associated charged particles relative 
to a 71° trigger with transverse momentum pj t > 7 GeV/c are shown for five intervals of associ- 
ated particle transverse momentum pj. In all cases, strong correlation peaks on flat backgrounds 
are clearly visible, indicating the di-jet structure which is contained in an interval A(p = ±60° 
about a direction towards and opposite the to trigger for all values of associated pj (> 0.3 GeV/c) 



shown. The width of the peaks about the trigger direction (Fig. |14|a), or opposite to the trigger 
(Fig. |j"4|b) indicates out-of-plane activity from the individual fragments of jets. The trigger bias 
was directly measured from these data by reconstructing the trigger jet from associated charged 
particles with pj > 0.3 Gev/c, within Atp = ±60° from the trigger particle, using the algorithm 
PTjet = Pt, + 1.5£prcos(A0), where the factor 1.5 corrects the measured charged particles for 
missing neutrals. The measurements of (ztrfg) = (pt,/ PTjet) as a function of pj t for 3 values of 
(Fig. |l5]-(left)) show a variation which is consistent with scaling as a function of xj, which was 



not expected [71, 12]. Another observation [24], not much emphasized at the ISR but relevant to 
recent observations at RHIC, is that the measured (ztrig) i s different for single particle inclusive 
triggers and pair triggers (Fig. 15 -(right)). A recent measurement by STAR at RHIC J72[] gives 
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Figure 14: a,b) Azimuthal distributions of charged particles of transverse momentum pr, with respect to 
a trigger 71° with pn > 7 GeV/c, for 5 intervals of pr'- a) (left-most panel) for A(j) = ±n/2 rad about the 
trigger particle, and b) (middle panel) for A<p = ±n/2 about % radians (i.e. directly opposite in azimuth) to 
the trigger. The trigger particle is restricted to |t}| < 0.4, while the associated charged particles are in the 
range | T7 1 < 0.7. c) (right panel) xe distributions (see text) corresponding to the data of the center panel. 



(z t ) « 0.78 ± 0.04 for inclusive p Tl =7.0 GeV/c at yfs = 200 GeV. 

Following the analysis of previous CERN-ISR experiments [f73|, 17], the away jet azimuthal 
angular distributions of Fig. |l4|b, which were thought to be unbiased, were analyzed in terms of the 
two variables: p out = pj sin(A0), the out-of-plane transverse momentum of a track; and %e, where: 



x E 



-PT-PTt 
\PTt\ 2 



-prcos(A0) 

PTt 



JET 



— (5.20) 

Ztng 




Ztrig c± pjt I Prjet is the fragmentation variable of the trigger jet, and z is the fragmentation variable 
of the away jet. Note that xe would equal the fragmenation fraction z of the away jet, for Zttig 1> 
if the trigger and away jets balanced transverse momentum. The xe distributions [[70L for the 



data of Fig. [T4p are shown in Fig. |14p and show the fragmentation behavior expected at the time, 
e~ 6z ~ e -6x E {ztng) _ jf w idth of the away distributions (Fig. [I4|b) corresponding to the out of 
plane activity were due entirely to jet fragmentation, then (|sin(A0)|) = (\JtA/pt) would decrease 
in direct proportion to l/pr, where yV is the component of jj in the azimuthal plane, since the 
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Figure 15: CCOR [|4|] measurement of (ztrig) as a function of pj f - (left) for inclusive 7T 's and (right) for a 



7T -pair trigger with m x o n o > 8 GeV/c 



jet fragmentation transverse momentum, jj, should be independent of pj. This is clearly not the 



case, as originally shown by the CCHK collaboration [17], which inspired Feynman, Field and Fox 



(FFF) [jllj] to introduce, kj, the transverse momentum of a parton in a nucleon. In this formulation, 
the net transverse momentum of an outgoing parton pair is \flkj, which is composed of two or- 
thogonal components, y/lkj^, out of the scattering plane, which makes the jets acoplanar, i.e. not 
back-to-back in azimuth, and \/2kj x , along the axis of the trigger jet, which makes the jets unequal 
in energy. Originally, kj was thought of as having an 'intrinsic' part from confinement, which 
would be constant as a function of x and Q 2 , and a part from NLO hard-gluon emission, which 
would vary with x and Q 2 , however now it is explained as 'resummation' to all orders of QCD [|75|]. 



FFF [Jllj, [76p gave the approximate formula to derive kj from the measurement of p out as a 
function of xe '■ 



(|Pout|} 2 =4[2(fe|) 2 + (|jVJ) 2 ] + (Ur, 



(5.21) 



CCOR [77] used this formula to derive (\kr t |) and (| jV |) as a function of pj t and y 7 ? from the data 
of Fig. |l4|b. This important result showed that is constant, independent of pr, and y/s, as 

expected for fragmentation, but that (l&rj) varies with both pr t and y/s, suggestive of a radiative, 
rather than an intrinsic origin for kj. The analysis was repeated, this year, by PHENIX for p-p 



collisions at y/s = 200 GeV [Tgj. 



5.3 Why 'everything' became 'almost everything' due to a new understanding of xe 
distributions 



The new measurement of (jr), (&r) and the xe distribution, this year [10], led to several 
complications and surprises, most notably that the shape of the xe distribution is not sensitive to 
the fragmentation function. The complications concern the fact that while the effect of (it) could 
be neglected at the ISR, where (z t ) ~ 1 due to the larger value of n, it had to be taken into account 



at RHIC, with the result that the already complicated formula (Eq. 5.21) for deriving kj became 
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even more complicated: 



(Zt(k T ,X h )} 

{xh(k T ,x h )) 



(5.22) 



where x/, (xh) is the ratio of the associated particle (parton) transverse momentum to the trigger 
particle (parton) transverse momentum: 



X/, = Xfj — XfryKj ,Xh) = — 

Pt, Pt, 



(5.23) 



Note that the hadronic variable Xh is measured on every event and that the partonic variable Xh is a 
function of both kj and Xh (Fig. 16 -(right)), as is the "trigger bias" (z, t ) [10]. Thus, the solution of 



Eq. 5.22 for y(k T ) is an iterative process. The results Jl0| ] are shown in Fig. 16 
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Figure 16: (left) y (k^} values in p-p collisions for associated charged particles with 1 .4 < pj B < 5 GeV/c 
as a function of pr, of a n° trigger (solid symbols) with statistical and systematic errors, measured by 
PHENIX |@. The CCOR masurement at V* = 62.4 GeV Q (open triangles) is also shown, (right) (z t ) 
and (xh) as a function of pj t used in the PHENIX measurement, shown with statistical and systematic errors. 



In order to evaluate (z t ) the fragmentation function must be known. Based on the longstanding 
belief that the away jet was unbiased, PHENIX attempted to derive the fragmentation function from 
the measured xe distribution. 

5.3.1 Fragmentation Formalism — two-particle correlations from a jet-pair 

First recall the joint probability for a fragment pion, with pj t = z t pj t , originating from a parton 



with Pt, (Eq. 5.1): 



d 2 o n (p Tt ,Zt) 



da,. 



xD«Jz t ) 



p Tf dp Tt dzt PT,dp Tl 

= f q (pT t )xD*(z t ) 



(5.24) 
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Here we make explicit that f q {pT,) represents the kj -smeared final-state scattered-parton invariant 
spectrum do q j pj t dpj t and D*(zt) represents the fragmentation function. Due to the kj smearing, 
the transverse momentum pj a of the away partem in the hard-scattered parton-pair is less than the 



transverse momentum of the trigger parton p Tt [JlOp . The probability that the partem with p Ta frag- 
ments to a particle with pj a = z a pT a in interval dz a is given by D™(z a )- Thus, the joint probability 
for a fragment pion with pj t = z.t Pt, , originating from a parton with pj t , and a fragment pion with 
Pr a = z a pT a , originating from the other parton in the hard-scattered pair with pj a is: 



d 3 o 7l (p T ,,Zt,Za) da q 



pT f dpT,dz t dza pT f dpT, 



xD*(z t )xD«Jz a ) 



(5.25) 



where 



Z a 



Pja_ 
PT a 



PTa 



ZtPTa 



XhPT, XhPT, 



and x/j = pjj Pt, (Eq. |5.23| ). Changing variables from pj t , z, t to pj t , Zt as above and similarly from 
z a to p Ta yields: 



d 3 e n 



1 



do,. 



-^(*)°J(??r) 



dp Tt dz. t dpT a XhPT, d(p Tr /Zt) ' '' XhPT, 

where for integrating over Zt or finding (z t ) for fixed pt,, PT a , the minimum value of Zt is zf 
2pT t /-y/s = xt, and the maximum value is: 



(5.26) 



max 

z t 



r. PT, 

Xh — 

PT a 



Xh 
Xh 



where Xh(PT„PT a ) is also a function of £7- (Eq. 5.23 ). Integrating over dzt in Eq. 5.26| gives the xe 
distribution in the collinear limit, where pt ci = xePt,, and it was thought [|K]] that a simply parame- 
terized fragmentation function could be extracted from a joint fit to the measured xe and inclusive 



Pt, distributions (Eq. 5.6). However, there were serious difficulties with convergence which took a 



while to sort out. Eventually, the xe distributions were calculated from Eq. 5.26 using LEP mea 



surements for quark and gluon fragmentation functions, with shocking results (see Fig. 17) — the 



xe distributions calculated with quark D q pa exp{— 8.2 • z) or gluon D n g « exp(— 11.4 -z) fragmen- 
tation functions do not differ significantly ! Clearly, the xe distributions are rather insensitive to the 
fragmentation function of the away jet in contradiction to the conventional wisdom dating from the 
early 1970's. 

The evidence of this explicit counter example led to an attempt to perform the integral of 
Eq. 5.26| analytically which straightforwardly confirmed that the shape of the xe distribution is not 
sensitive to the shape of the fragmentation function. However, it was found that xe distribution 
is sensitive to Xh, the ratio of the transverse momentum of the away-side jet (pT a ) to that of the 
trigger-side jet (Pt,)- This can be put to use in A+A collisions to measure the relative energy loss 
of the two jets from a hard-scattering which escape from the medium. 

5.3.2 Analytical formula for the xe distribution 



With a substitution of a power-law parton pj t spectrum (Eq. |5.5| ) and an exponential fragmen- 
tation function (Eq. 5.7 ), as in section |5.1.1| , the integral of Eq. 5.26 over z.t becomes: 



d 2 o n 
dpT,dp Ta 



B 2 A 

Xh P\ J *t, 



P Ta , Ti-1 r 

dz t z t exp[ 



-bztiX + P^-)] 

XhPT, 



(5.27) 
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Figure 17: (left) xe distributions from PHENIX [Il0|] in p-p collisions at ^fs = 200 GeV for several values 



of pt, ■ The solid and dashed lines represent calculations of the distribution from the integral of Eq. 5.26 for 
quark (solid lines) and gluon (dashed lines) fragmentation functions based on exponential fits to the LEP 
measurements [J78L \l% shown on the right panel. 



This is again an incomplete gamma function, if Xh is taken to be constant as a function of Zt for 
fixed p Tl , pr a - 

d 2 o n B 1 A 1 



dpT,dp Ta Xh p\ b ' n 



T(n, b'xj ) — r(«, b'xf, — — ) 



where b' is given by: 



b'=b(l + ^ 



XhPT, 



(5.28) 



(5.29) 



The conditional probability of the pj a distribution for a given pj t is the ratio of the joint prob- 
ability Eq. 5.28 to the inclusive probability Eq. Q3] , or 



dpi 



B 



Y{n,b'x Tt )-T{n,b'x h ^- a ) 



PTl 



bp Tt x h (1 + ^ Y [T{n - 1 , bx Tt ) - r(n - 1 , b)] 



(5.30) 



and this answer is exact for the case of constant x/ u with no assumptions other than a power law for 
the parton pj t distribution and an exponential fragmentation function. In the collinear limit, where, 

PT a =X E PT,- 



dPn 

dx£ 



1 B 



1 



r(n,b'x Tt )-r(n,b'f E ) 



PT, 



Xh b{\ + fY \T(n-l,bxT t )-r(n-\,b)] 



(5.31) 



With the same approximation for the incomplete gamma functions used previously (Eq. |5.17| ), 
namely taking the upper limit of the integral (Eq. 5.27 ) to infinity and the lower limit to zero, the 
ratio of incomplete gamma functions in Eq. 5.31 becomes equal to n — 1 and the xe distribution 
takes on a very simple and very interesting form: 



dX£ 



(m) (n-i; 



1 



1 



Pt, 



x h {i + x ±Y 



(5.32) 
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where the only dependence on the fragmentation function is in the mean multiplicity of charged 



particles in the jet (m) mB/b^b. The dominant term in Eq. |5.32| is the Hagedorn function 1/(1 + 



xe /xti) n so that Eq. 5.32| exhibits x^-scaling in the variable xe /xt,. The shape of the xe distribution 



is given by the power n of the partonic and inclusive single particle transverse momentum spectra 
and does not depend on the exponential slope of the fragmentation function. However, the integral 
of the xe distribution (from zero to infinity) is equal to (m), the mean multiplicity of the unbiased 
away-jet. 

The reason that the xe distribution is not very sensitive to the fragmentation function is that the 



integral over z.t for fixed pj t and pj a (Eqs. 5.26, 5.27) is actually an integral over the jet transverse 



momentum pj r However since both the trigger and away jets are always roughly equal and opposite 
in transverse momentum, integrating over pj t simultaneously integrates over pj a , and thus also 
integrates over the away jet fragmentation function. This can be seen directly by the presence of Zt 



in both the same and away fragmentation functions in Eq. |5.26| , so that the integral over z.t integrates 
over both fragmentation functions simultaneously. 

5.3.3 Why did we believe that the xe distribution measured the fragmentation function? 

The seminal paper of Feynman, Field and Fox (FFF) [ [TT| ] was especially influential in form- 
ing the belief that the xe distribution measured the fragmentation function. To cite directly from 



Ref. [11], p 25, "There is a simple relationship between experiments done with single-particle trig- 
gers and those performed with jet triggers. The only difference in the opposite side correlation 
is due to the fact that the 'quark', from which a single-particle trigger came, always has a higher 
p± than the trigger (by factor l/zt). The away-side correlations for a single-particle trigger at p± 
should be roughly the same as the away side correlations for a jet trigger at /?^(jet)=/?^(single 
particle)/ (&)". This point is reinforced in the conclusions (p 59), "2. The distribution of away- 
side hadrons from a jet trigger should be the same as that from a single particle trigger except 
for a correction due to (z t ) (see Fig.23)" [which is shown as Fig. |TJ below]. Another interesting 
point is, "8. Because the quarks scatter elastically (no quantum number exchange - except perhaps 
color), the away-side distribution of hadrons in pp collisions should be essentially independent of 
the quantum numbers of the trigger hadron." — i.e. the jets fragment independently. Note that in 



FFF the notation is a + b — > c + d (as in Eq. |2.1| ) where a,b,c,d are called 'quarks' , so FFF call Zt, 

z c - 

This belief was thought to have been verified by measurements at the CERN-ISR which 
showed (Fig. [T^-(left)) that jet fragmentation functions in v-p, e + e~ and p-p reactions (CCOR 
Fig. 14: [ f70| , [74]]) are the same, with the same dependence of the exponential slope b on s (Fig. 19- 
(right) @]. 

5.3.4 A very interesting formula 



Equation Eq. 5.32 (repeated below in a slightly different format) is very interesting. 



dy 



7a(m)(n— 1)-. — — r- where y=-r- ■ (5.33) 



It relates the ratio of the transverse momenta of the away and trigger particles, pj a j pj t = %h « xe, 
which is measured, to the ratio of the transverse momenta of the away to the trigger jet, Pr a /PT,, 
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Fig. 23, Comparison of the tt-Q and jet tnyijei iiway-sidu distribution of charged hadron&in pp 
collisions at W * 19.4 CeV, 6 s = 90°, wndpi {trigger) ■ 4.0 GeVfc from the quark-quark scat- 
ter ing mode!. The upper figure shows the single-particle trigger results pint ted versus 
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Figure 18: Figure 23 from FFF [[j]} 



: 0.2 
0.05 

















i i i i 



0.2 03 0.4 O.S D.S 07 

Figure 2! Jl-I l'rii:'.nii'-irn-iiir: fi;n:'liivi<i ri'fi.K :i : .d in Ji::erei;L pruCrKvis : ■.'.-> iniEraetioTisfcpsj-l 
trinn^es, Van der Welde 1979}; e + e~ an nihil nitons {solid tine, Hanson -el ai 1 97 j]: and pp 
c&ltistonsl/uJIdrdfiCS,^ < tGsV!c, open circles CS,p^. > 6 GeV / e, /utf squares CCOJt.jJt 
> 5 GeV/c, of^ii^itfrres CCGR./h > 7 GcV/£j. 



km 



ohigh jet 
• e + e" (Spear) 



J_ 



2 4 6 pjGeWfc 

Figure 1 9 Tht slopes b obtained from CXiXjnential Gls to the jet fragmentation function in lire 
interval 6,!<?f 0,8 in s*e" annihilation (full rirclm) and LPTH data ot the BS 
Collaboration {open circles). 



Figure 19: (left) Jet fragmentation functions from v-p, e + e and p-p reactions, (right) ^-slopes from this 
data, where 'LPTH' is an acronym for Large pr Hadron production [E7I]. 
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which can thus be deduced. Although derived for p-p collisions, Eq. 5.32 (5.33) should work 
just as well in A+A collisions since the only assumptions are independent fragmentation of the 
trigger and away-jets with the same exponential fragmentation function and a power-law parton 
Pt, distribution. The only other (and weakest) assumption is that x/, is constant for fixed pj t as 
a function of x#. Thus in A+A collisions, Eq. 5.32 for the xe distribution provides a method of 
measuring the ratio x/, = Pt c J Pt, and hence the relative energy loss of the away to the same side jet 
assuming that both jets fragment outside the medium with the same fragmentation function as in 
p-p collisions. 



5.3.5 Test of determination of Xh from the xe distribution in p-p collisions 




Figure 20: (left) Eq. 5.32 for n=8.1 divided by (m). The integral should be equal to 1. Curves are for 
i/, =1.0 (red), 0.8, 0.6, 0.4, 0.2 (blue), with intercept=7.1/i /l at x E = 0. (right)-PHENIX x E distribution for 
5 < p T , < 6 GeV/c (Fig. ^(left)) with Eq. |3| for x h = 0.8 (black); STAR [|72| x E distribution (Fig ||) 



with Eq. 5.32 for i/, = 1.0 (blue). Curves are from left panel multiplied by a factor of 1.63 to agree with 
PHENIX data. STAR data have been multiplied by a factor of 0.6 to agree with this normalization. 



A plot of Eq. |5.32 is shown in Fig. |2^-(left) for n = 8.1 for various values of X},. Clearly, the 
smaller the value of x/,, the steeper is the xe distribution. However, all the curves in Fig. ^Oj-(left) 
are related by a simple scale transformation of Eq. 5.33 : y — > xe = Xhy- In general, the values 



of n and x/, should be able to be determined from a simultaneous fit of the inclusive pj t spectrum 



(Eq. [331 or p.l5[ ) and the xe distribution (Eq. |5.31| ). On the other hand, when the value of n is 
well determined from the inclusive pj t spectrum (e.g. see Fig. Eq. |5.33 can just be scaled to 
fit the measured xe distribution. This was done for the PHENIX p-p data JT0| ] of Fig. |l7|-(left) 
(5.0 < pt, < 6.0 GeV/c). The value of x/ 2 = 0.8 which gave the best "eyeball" agreement with the 



data (Fig. 20-(right)) agrees with the value of x/, determined independently from the ^-smearing 
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analysis (Fig. 16) [1C]. Similarly, for the STAR p-p data [72] (Fig. 25) excellent agreement is 
found withx/, = 1.0, for xe > 0.2. Thus, the method works for p-p collisions. Note that the evident 



deviation of the STAR data from Eq. |5.32| for xe < 0.2 may be a limitation due to the simple 
approximations or may be the result of the absence of corrections for decay-in-flight for the low 
pj a particles (presumably in the measurement. Also, as indicated in the caption of Fig. |2(| 



typically, Eq. 5.32 must be normalized in order to agree with the data. 



5.4 2-particle correlation measurements in Au+Au collisions 

One of the first and still most striking measurements of two-particle correlations in Au+Au 
collisions was presented by STAR at the Quark Matter 2002 conference [|80|]. In Fig. 21, the 
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Figure 21: STAR [ |80| ] conditional probability j^—-^ °f associated non-identified charged particles with 
Pt„ in the range 2GeV/c < pj a < pj t per trigger particle with pj t between 6 and 8 GeV/c, all in the range 
< | At] | < 1.4 for central Au+Au collisions at ^/snn = 200 GeV as a function of the azimuthal angle 
difference A0 between the trigger and associated particles (solid circles). Non-jet background modulated by 
elliptic flow v? is shown as the blue curve. The sum of the measured p-p correlation plus the flow is shown 
as the red histogram. 

conditional probability — given a trigger particle with pj t between 6 and 8 GeV/c — of detecting an 
associated particle with pj a in the range 2GeV/c < pj a < pj t is shown for central Au+Au collisions 
at y/sNN = 200 GeV as a function of the azimuthal angle difference of the two particles, A(j>. The 2- 
particle correlation function expressed as the conditional probability is the sum of the background 
of particles randomly associated to the trigger, which is modulated by the common hydrodynamic 
flow (represented by V2(pr)), plus the jet correlation function which was presumed to be the same 



as that measured in p-p collisions [81]: 

'(A0)=Cf (A0)+*(l + 2v2(pr,)v2(p r Jcos(2A0)) 



/-iAkAh / 
c 2 



(5.34) 



The trigger-side correlation peak in central Au+Au collisions appears to be the same as that mea- 
sured in p-p collisions (corrected by the small flow effect) but the away side jet correlation in 
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Au+Au appears to have vanished, the data seem to be saturated by the small flow effect. This obser- 
vation appears consistent with a large energy loss in the medium, or a medium that is opaque to the 



propagation of high momentum partons, as originally indicated by the suppression observed [ |82| ] 
in single particle inclusive measurements for pj t > 3 GeV/c (recall Fig. |5|). 

Although the apparent vanishing of the away jet in central Au+Au collisions is fantastic from 
a public relations perspective, it is misleading from a scientific viewpoint as it suggests that the 
away-jet was totally absorbed by the opaque medium. Later work presented by STAR at Quark 
Matter 2004 [I 



\ |72(] with 4 < p T> < 6 GeV/c and 0.15 < pr a < 4 GeV/c showed that the away jet 
didn't disappear, it just lost energy and the away-side correlation peak became much wider than in 



p-p collisions (Fig. 22). 
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Figure 22: STAR conditional probability 2-particle correlation function with flow-modulated background 
subtracted: (left) Measurements in d+Au [0], p-p and Au+Au central |j8^| collisions at ^/snn = 200 GeV 
with 4 < pj, < 6 GeV/c and 2GeV/c < pj a < pj t \ (right) STAR data with same trigger pj t but with 0.15 < 
p Ta < 4 GeV/c & 



Still later work presented by STAR at Quark Matter 2005 [86, 87], this past year, showed that 
an away-side jet correlation peak with the same width as in p-p collisions re-appeared when pj t 
was raised to the range 8 < pj t < 15 GeV/c, with pj a > 3 GeV/c (Fig. Q). Clearly, the study of jet 
phenomena by two-particle correlations in Au+Au collisions is much more complicated than the 
same subject in p-p collisions and one can expect a long learning curve. 

However, even at this early stage, there is one definitive result from 2-particle jet correlations 
(Fig. |4j) p9|], in the sense that it casts serious doubt on the explanation of the 'baryon anomaly' 
(recall Fig. 0) by coalescence models [ 46 , ^7 , fjq ]. Fig. |24"|-(left) shows the conditional probabil- 
ity 2-particle azimuthal correlation functions, with the integrated associated particle yields/trigger 
shown in Fig. 24 -(right), for p-p, d+Au and AuAu collisions in which the trigger is either an iden- 
tified meson or baryon in the range 2.5 < pj t < 4.0 GeV/c and the associated particles, in the 
range 1.7 < pr a < 2.5 GeV/c, are not identified. The yield of associated particles/per trigger on the 
near side, from the same jet as the trigger hadron, is the same for meson and baryon triggers as a 
function of centrality, except perhaps in the most central bin; and the same effect is seen for the 
away-side yields. The red-dashed curve indicates the expected trigger-side conditional yield if all 
the anomalous protons in Au+Au collisions were produced by coalescence. This shows that meson 
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Figure 23: STAR conditional probability 2-particle azimuthal correlation histograms for charged hadron 
triggers with 8 < pj, < 15 GeV/c, in minimum-bias d+Au, 20-40% Au+Au and 0-5% Au+Au collisions at 
v/*aw=200 GeV. pr a increases from top to bottom as indicated pft p7\\. 
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Figure 24: (left) PHENIX [ 49 1 conditional probability 2-particle azimuthal correlation distributions ^J— 



for triggers by identified mesons and baryons with 2.5 < pr, < 4 GeV/c and associated non-identified charged 
hadrons with 1 .7 < pj a < 2.5 GeV/c in collisions at ^/snn = 200 GeV: Au+Au, as a function of centrality as 
indicated, d+Au minimum bias and p-p (non-identified triggers), (right) Conditional yield per trigger meson 
(circles), baryon (squares) from this data, integrated within A(j) = ±0.94 radian of the trigger (Near Side) or 
the opposite azimuthal angle (Away Side), for Au+Au (full), d+Au (open) collisions at ^/s^n = 200) GeV. 
Shaded boxes indicate centrality dependent systematic errors. An overall systematic error which moves all 
the points by 12% is not shown, p-p data are shown for non-identified charged hadron triggers. 
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and baryons at intermediate pj are produced by hard-processes with the same di-jet structure, and 
not by soft coalescence. 



5.4.1 Jet energy loss or jet absorption? 

In section 5.3.4] , I asserted that the ratio of the transverse momenta of the away-jet to the 
trigger-jet, x\ r = pT a /pT t , and hence the relative energy loss of the away to the same side jets in 
both p-p and A+A collisions could be determined from measurements of the xe distribution using 
Eq. |5.32| . A test of this method, which worked for p-p collisions, was presented in section 5.3.5 
Now I apply the method to Au+Au collisions. 



Fig |25"l-(left) shows the STAR measurement [72] of the pj a distribution, given pj,, from the 
data shown in Fig. |22|-(right). In Fig. 25 -(right), these measurements are plotted as an xe ~ Pt„/pt, 




* STAR* PRL95 p+p 

□ STAR* PRL95 AuAu 40-80 

* STAR* PRL95 AuAu 00-05 




Figure 25: (left) STAR measurement [ f72| ] of transverse momentum (p±) distribution of associated charged 
hadrons for a trigger charged hadron with 4 < pr, < 6 GeV/c for pp, Au+Au peripheral(80-40%), Au+Au 
central (top 5%) collisions at </Smv = 200 GeV. a) near-side, b) away-side c,d) /AA=ratio of AA to pp p± 
distributions for c) near side, d) away side, (right) data from (b) plotted as dP/dxE compared to Eq. p".32 
with x.h = 1 for p-p, xi, = 0.75 for Au+Au peripheral, and x/, = 0.48 for Au+Au central. The normalization 
is the same as in Fig. ^-(left) for p-p collisions for both the data and the curve. The Au+Au data have been 
multiplied by the same factor of 0.6 to maintain the relative normalization as published, but the curves are 
normalized and i/,-scaled by eye to agree with the measurements in the range 0.2 < xe < 0.8. 



distribution and shown together with Eqs. |5.32| , |5.33| , with n = 8.1, scaled to match the data, which 
are beautifully consistent with no relative energy loss of the two jets in p-p collisions as noted above 



(recall Fig. |20j-(right)). By contrast, in Au+Au collisions, agreement with the data is obtained with 
a ratio of away/trigger jet momenta of 0.75 in peripheral (40-80%) and 0.48 in central (0-5%) 
collisions. This indicates a clear relative energy loss of the away jet compared to the trigger jet, 
which increases with increasing centrality. However, the trigger jets in Au+Au are surface biased 
by the falling power-law pj spectrum, an effect analogous to 'trigger bias' — the jets which give 
trigger particles of a given pr t are more likely to be produced near the surface and lose little energy, 
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with pr t close to pj t , than to have been produced deeper in the medium with a larger pj t 3> pr, 
and then have lost significant energy in getting to the surface. Because of the trigger-jet surface 
bias, the away-jets must traverse the entire medium in order to be observed (except for the unlikely 
cases when the jet-pair is tangential to the medium). Hence, the decrease in x% from 1.0 in p- 
p collisions to 0.75 in Au+Au peripheral (40-80% ) collisions to 0.48 in Au+Au central (0-5%) 
collisions indicates that the energy loss of the away-jet increases with distance traversed in the 
medium. 



I then tried to analyze the higher pj t STAR away-jet measurement [gg, |87j] by the same 



method [88]. First, I plotted the data from both STAR measurements as xe distributions on the 
same scale (Fig. l26|). The measurements appear to disagree, both in normalization and shape, so I 
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Figure 26: STAR data [|72|] from Fig. ^ for d+Au, Au+Au (peripheral) and Au+Au central (0-5%) collisions 
at y?SNN = 200 GeV (filled points) compared to higher pj t STAR data [p7l p8|| (open points). 



tried normalizing the higher pj t measurement [ p7[ ] to agree with the lower pj t measurement Q72], 



which would be correct if xe scaling were valid in Au+Au collisions (see Fig. 27). The results are 
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Figure 27: Data from Fig. with the higher pr t data normalized to agree with the lower pj t data for 
xe < 0.45, together with the curves from Fig. [25] which fit the lower pj t data. 



quite interesting. It appears that the points at lower xe for the higher pj t measurement are consistent 
with the shape of the lower pj t data for xe < 0.45, with a dramatic break and a flattening of the slope 
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for xe > 0.5. This could be suggestive of a two-component distribution where some jets, which 
pass through the medium, lose energy, while other jets, such as those emitted tangentially, punch 
through without any energy loss. However it is difficult to understand why the punch-through of 
tangential jets would depend on the trigger pj r The comparison of the two STAR measurements 
and the possibility of a dramatic break in the xe distribution would be greatly clarified if a few 
lower xe points could be obtained for the higher pj t data. 
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Figure 28: (left) STAR higher pr, data [p7j|88|] as plotted in Fig. |2q together with curves of Eq. p. 32| adjusted 
by eye to best represent the data, (right) STAR higher and lower pj t data from Fig. for Au+Au central 
collisions together with PHENIX p-p data [|o[ from Fig. 20 (right). 



It is also possible to compare Eq. |5.32| to the STAR higher pj t data without reference to the 
lower pj t data (see Fig. ^8|-(left)). Here, another troubling effect is revealed. The best values of £/, 
are 1.30 for the d+Au data, 1.20 for the Au+Au peripheral (20-40%) data and 0.85 for the Au+Au 
central (0-5%) data. Thus, to within the error of the simplistic "eyeball" scaling, the away-jet in 
Au+Au central collisions with higher pj t , pj a also seems to lose about half it's energy relative to 
d+Au, consistent with the lower pj t measurement. However the xe slope for the higher pj t data 
is much flatter than other measurements in p-p and d+Au collisions in the same pj t range (see 
Fig. 



-(right) and Fig |29|) as reflected in the anomalous value of x/, = 1.30. In any model of jets 
with kj smearing, £/, must be < 1 as indicated by the other STAR and PHENIX data at RHIC 
(Figs. ^(fl-(right), 28 -(right), 29 -(right)). This clearly warrants further investigation. 



5.5 Possible new effects revealed by correlation measurements in Au+Au collisions at RHIC 
5.5.1 The ridge 

Due to the large acceptance of the STAR detector, a near-side correlation in pseudo-rapidity 
(rj) covering the full STAR rj acceptance was detected in addition to the flow modulated back- 



ground and the near-side jet correlation pQ]. As indicated schematically in Fig. 30, the width of 



the ridge in the direction is comparable to the near-side jet correlation and must be taken into 



account in extracting the near-side jet yields. See Ref. [90] for details 
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Figure 29: PHENIX measurement | |8^ ] of conditional yield as a function of xe for near- and far-side n ± — 
correlations for several values of trigger pj t (indicated) in minimum bias d+Au collisions at ^snn = 200 
GeV. 
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Figure 30: (left) STAR |Q extraction of the near-side jet yields for 3 < pr, < 4 GeV/c and p Ta > 2 GeV/c: 
Jet+Ridge, when the whole AT7 acceptance is included; Jet(A0), Jet(ATj), when account is taken of the ridge 
by two different methods, (right) Schematic drawing of the extent of the same-side ridge and jet correlation 
in the A77 and A0 directions. 



5.5.2 Wide jets and/or Mach Cones-2 particle correlations 



Measurements of non-identified charged hadron correlations by both PHENIX [ |91| ] and STAR 
in the "intermediate pr" region (where the 'baryon anomaly' is found) are shown in Fig. [H]. For 
both PHENIX and STAR the trigger and associated particles have 2.5 < pj t < 4.0 GeV/c and 
1.0 < pr a < 2.5 GeV/c but the PHENIX range in pseudorapidity is |t]| < 0.35 for both particles 
while for STAR |t]| < 0.7 for the trigger particle and 1 77 1 < 1 for the associated particles. Although 
disagreeing in absolute value, presumably due to the different v\ acceptances, the PHENIX and 
STAR measurements both exhibit a striking widening of the away-side correlation in going from 
peripheral to central collisions, with a strong hint of a local minimum (dip) developing at Atp = % 
for centralities less than ~ 60%. The existence of these local minima per se is not significant once 
the systematic errors on V2 are taken into account but it is clear that all the away-side distributions 
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Figure 31: Conditional yields of associated particles with 1.0 < pj a < 2.5 GeV/c per trigger particle with 
2.5 < pr, < 4.0 GeV/c as a function of A(j) for various centralities in Au+Au collisions at JJ^n =200 GeV 
after subtraction of the flow-modulated non-jet background: (left) PHENIX [^TJ data where both 1 and 2 
unit bands are shown for the (rms) systematic uncertainty in V2\ (right) STAR [|2|] data with systematic 
uncertainties indicated as histograms. 



in the more central samples for both PHENIX and STAR have a very different shape than in the 
most peripheral sample and all seem to exhibit a dip at A(j) = %. 

5.5.3 Deflected jets and/or Mach Cones-3 particle correlations 

There are numerous explanations for the possibly two-peaked structure, roughly 1 radian away 



from n, in the away-side distributions shown in Fig. |31], of which two are commonly discussed: 
a 'Mach cone' [[)3|] due to the away parton exceeding the speed of sound in the medium and 
causing the QCD equivalent of a sonic-boom; or deflected jets, due to the strong interaction with the 
medium which, e.g. for mid peripheral collisions where the overlap region has a large eccentricity, 
might prevent directly back-to-back jets from penetrating through the medium (see Fig. f32|-(left)). 
Both STAR and PHENIX try to distinguish a Mach cone from deflected jets using 3-particle 



correlations. In Fig. 32 STAR [94] studies the correlation of a trigger particle with 3 < pj t < 4 
GeV/c to 2 associated particles with 1 < pj a < 2 GeV/c by making a plot of A(j>\ versus Afa, 
the difference in azimuth of each associated particle with the trigger particle. In Fig. |2|-(center) a 
schematic of the expected results are shown on top for the case of deflected jets for which A<pi = A02 
when both associated particles are on the away-side and where A(p when one or both of the 
associated particles are on the trigger-side. The diagonal elongation near (n, 7i) is consistent with 



kj smearing, since the typical fragmentation transverse momentum J at 0.6 GeV/c is much 
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Figure 32: STAR method of 3 particle-correlations ||94j]: (left) Schematic drawing of deflected jets (top) 
and Mach cone (bottom) opposite to a trigger; (center) A schematic plot of the azimuthal angle differences 
A0i versus A02 of each of the two associated particles with respect to the trigger for deflected jets (top) 
and a Mach cone (bottom); (right) STAR measurement of A0i versus A02 in central Au+Au collisions at 
V^viv = 200 GeV. 



less than J (kj) 2.7 GeV/c [10], so that A0i « A0 2 / ft- For away-side particles which form 
a cone roughly around the direction opposite to the trigger, there are off-diagonal as well as on- 
diagonal A0i = A02 correlations (bottom). This is not obviously the best projection to understand 
this problem as illustrated by the measurement in Fig. |2|-(right) which is difficult to understand 
but does appear to show off-diagonal activity. 

PHENIX [^5|] defines a coordinate system for the correlation of two associated particles (1.0 < 
Pt ci < 2.5 GeV/c) to a trigger particle (2.5 < pr, < 4.0 GeV/c) in which a conical correlation would 
be directly visible (Fig f33|b). The angle 0* represents % — 0', where d' would be the half-angle 
of a cone centered opposite in azimuth to the trigger pj t \ and the variable A(j)* for this analysis 
represents the 'azimuthal' angle around the cone. The data are displayed as a polar plot of A(p* as 
a function of 6* (Fig. |33|a,c,d). For a Mach cone, there should be a '5-function' at a fixed half- 
angle 6* smeared by kj, and a uniform distribution in A<p* (Fig. 33c), while a deflected jet would 
show correlations that are close in both 6* and A(p* and would favor orientations of A(p* in the 
77 direction, since the same-side and away-side jets are relatively uncorrelated in pseudorapidity 
(Fig. 33 a). The measurement shown in Fig. 33d seems to exhibit both types of activity. 



6. Conclusion 

Much has been learned, both in the 1970's and recently at RHIC, by the study of jets and 
hard-scattering via single particle, two-particle and 3-particle measurements in p-p and A+A 
collisions. Clearly, for measurements in A+A collisions at RHIC, we are still at the early stages of 
a long and interesting learning process. 
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Figure 33: PHENIX method of 3-particle correlations [p5[], b) Schematic illustration of the coordinate 
system used for 3-particle correlations. The high pj, trigger particle serves as the near-side jet axis and 0* 
represents n — 0\ where 0' would be the half-angle of a cone centered opposite in azimuth to the trigger 
/%; the variable A<j>* represents the 'azimuthal' angle around the cone, a) Simulated 3-particle correlation 
for a deflected away-side jet in the PHENIX acceptance, c) Same as (a) but for a conical away-side jet. 
d) Measurement of 0*, A(j)* 3-particle correlation surface for charged hadrons in central (0-5%) Au+Au 
collisions within the PHENIX acceptance. 
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